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ABSTRACT 

In this Letter we propose a simple thermal comptonization model to account for the observed 
properties of the phase lags associated to the "plateau" intervals of GRS 1915+105. By invoking 
a temperature stratification in a corona and assuming that the optical depth of the comptonizing 
region increases as the disk inner radius moves inward, we are able to reproduce both the 
observed colors and time lags in the continuum. 



Subject headings: accretion, accretion disks — 
individual (GRS 1915+105) 

1. INTRODUCTION 

The galactic microquasar GRS 1915+105, discov- 
ered i n 199 2 by GRANAT (Castro-Tirado, Brandt & 
Lund 1992| ), is a transient source with an extremely 
wide variety of variability modes (Greiner et al. [L996| ; 
Belloni et al. 1997; Muno, Morgan & Remillard 
1999|) . During the period 1996-1997 it has been ex- 
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val of frequencies, ranging from about 0.001 to 67 Hz 



tensively and repeatedly observed by the Rossi X- 
ray Timing Explorer (RXTE), alternating phases of 
dramatic variability to phases of remarkably regular 
behavior (Belloni et al. |2000|) . 

GRS 1915+105 was the first galactic source to show 
superluminal radio expansion (Mirabel & Rodriguez 
1994 ), usually interpreted in terms of emission from 
gas expanding at relativistic velocities. Measure- 
ments of this superluminal motion allow to estimate 
the source distance, which is in the range 6-12.5 kpc 
(Rodriguez & Mirabel [L999| ; Fender et al. |1999Q . Be- 
cause of the high extinction along the line of sight, 
only a faint infrared counterpart has been discovered 



(Morgan, Remillard & Greiner 1997; Chen, Swank & 
Taam 1997 ). In particular, in some observations a 
QPO with variable centroid frequency {vqpo ~ 0.5- 
10 Hz), whose value and root mean square amplitude 
correlate with the thermal flux, is detected (Mark- 
wardt, Swank & Taam 1999] ; Trudolyubov, Churazov 
& Gilfanov p99|) . 

The X-ray spectrum of GRS 1915+105 shows all 
the distinctive properties of other Galactic Black Hole 
Candidates (BHCs). A power-law tail with variable 
luminosity and slope is present at high energies, while 
a thermal component often dominates at low ener- 
gies (Belloni et al. 1997 ; Muno, Morgan & Remillard 
1999j ). A detailed spectral classification of the differ- 



(Mirabel et al. 1994). As a consequence, neither 
the mass function nor the orbital period are known. 
The source is believed to host a black hole because 
it shows properties (variability and superluminal mo- 
tion) very similar to the other known galactic micro- 
quasar, GRO 1655-40, for which a precise dynami- 
cal measurement yields M ~ 7M (Orosz &: Bailyn 



1997| ) 



Quasi Periodic Oscillations (QPOs) have been ob- 
served in GRS 1915+105 in an incredibly large inter- 



ent states of GRS 1915+105 has been recently pre- 
sented by Belloni et al. (2000). According to these 
authors, the variability of the X-ray emission is as- 
sociated with the alternation of three basic states, 
called A, B and C (see also Markwardt, Swank & 
Taam [TM|; Muno, Morgan & Remillard ^999| ). Dur- 
ing states A and B, corresponding to relatively high 
and soft flux, a strong thermal disk component dom- 
inates the emission. In state C, the 3-25 keV flux 
is characterized by lower luminosity and harder X- 
ray colors, being dominated by a hard (r = 1.5-2.5) 
power law with much softer (or even unobservable) 
thermal component. State C can last uninterrupt- 
edly, for a long time (weeks to months, the so called 
plateau intervals), in which case no soft component 



2 



PHASE LAGS IN GRS 1915+105 



is observed at all. State C is always associated with 
the 0.5-10 Hz QPO and a flat-top component in the 
Power Density Spectrum (Markwart, Swank & Taam 
1999| ; Trudolyubov, Churazov & Gilfanov 199S; Rcig 
et al. pOOOp . 

The analysis of a set of RXTE observations of GRS 
1915+105 in this state showed that phase lags are 
present, both in the continuum and the QPO (Reig 
et al. pOOOp . The sign of the lag is correlated with 
QPO frequency, count rate and X-ray colors. The lag 
is positive (photons in the 5-13 keV band arrive later 
than those in the 2-5 keV band) when the spectrum 
is hard and the centroid of the 0.5-10 Hz QPO is be- 
low 2 Hz; it then becomes negative (soft photons are 
delayed) as the spectrum softens and the QPO cen- 
troid frequency increases above 2 Hz. The smooth 
correlations of the lags with colors suggests that a 
unique mechanism is driving the observed transition 
from positive to negative lags. 

In this letter we propose a simple comptonization 
model to account for the observed properties of GRS 
1915+105 in state C, which is outlined in § ||. In 
§ H we report the main results obtained and compare 
them with observations. Finally, § [| contains a short 
discussion of the implications of our results. 

2. THE MODEL 

No definite model for the accretion flow in GRS 
1915+105 has been presented as yet. Here we con- 
sider a highly idealized scenario in which a Shakura- 
Sunyaev disk (Shakura & Sunyaev 1973] ) coexists 
with a hotter component. As suggested by Bel- 
loni et al. ( 1997] ), the main cause for the large- 
scale time variability of this source is the onset of a 
thermal-viscous instability which blows off the inner, 
radiation-pressure supported region of the disk. In 
particular, according to this scenario, the hard state 
corresponds to the replenishment by the disk of the 
previously evacuated region on a viscous timescale. 
Further support to this picture has been provid ed by 
numerical simulations (Szuszkiewicz &z Miller [1997 ; 
1998| ). 

Fitting the observed spectra of BHCs in the hard 
state requires the presence, in addition to a standard 
disk, of a hot, less dense phase, formed e.g. by the 
debris of the puffed-up inner region or by the evap- 
oration of the disk itself. The upscattering of soft 
disk photons by a hot cloud at T « 10 keV can ac- 
count also for the positive lags of the continuum. In 
fact, Compton cloud models at constant temperature 
have been already proposed in connection with other 
sources which show (positive) time delays (Miyamoto 
et al. 1991; Vaughan et al. 1994; see also Pouta- 
nen 2000 for a review). Compton models suffer from 



the problem that soft photons cool down the cloud 



to the temperature of the impinging radiation, un- 
less an efficient heating mechanism is at work. How- 
ever, if cooling of ions is inefficient, like in advection- 
dominated flows, protons decouple from electrons and 
stay hot at a fraction of their virial temperature, 
T p {r) ~ rjGMrrip/r. The electron thermal balance 
is fixed by the competition between ion heating and 
Compton cooling. A simple estimate of the equilib- 
rium electron temperature gives 

T(r) « 10 {c 2 r m /GM) 2 '\l - r in /r)- 2 / 5 keV (1) 

for a scattering depth ^ 1, rj = 0.01, a nd L s o f t oc 
(1/rtn — 1/r). As expected (see Poutanen 2000|) , the 
cloud temperature diminishes as the disk rim moves 
in and L so f t increases, but the dependence is weak. 
In an ion-heated corona the temperature is larger at 
smaller radii because L so ft decreases with r so that 
less photons are available to cool the inner part of the 
cloud. 

A non-homogeneous Compton cloud model may in- 
deed explain the positive/negative lags and the de- 
crease in the hard color observed in the state C of 
GRS 1915+105, provided that the corona becomes 
denser as the disk extends inwards. The density may 
increase in time if a cloud of fixed mass contracts as 
the result of angular momentum losses caused by vis- 
cous and/or radiation stresses. In order to illustrate 
the basic idea of the model, we sketch the corona 
as formed by two regions (referred to as "warm" and 
"hot"), each at uniform density and temperature, T\y 
and Tu , with Tjj > Ty/ (see Figure [j]) . The scattering 
depths are T\y and th, th > tw- When the density 
of the corona is low, soft photons get preferentially 
scattered in the inner part of the corona which has 
T H ^ 1 (while tw ^1), gaining energy and produc- 
ing a positive lag. As density increases, the Comp- 
ton y-parameter of the corona exceeds unity below 
a certain radius r sc . Photons emitted by the inner 
part of the disk undergo repeated scatterings with 
the hot electrons and start to fill the Wien peak at 
Th- The spectrum emitted by the disk plus the hot, 
dense part of the corona is therefore the superposi- 
tion of a multitemperature blackbody and a Bose- 
Einstein distribution at Th] the latter has the same 
number of photons originally emitted by the disk in 
the region rj n <J r <^ r sc . These photons have still 
to traverse the warm part of the corona before escap- 
ing to infinity. In doing so, they are downscattered 
and this produces the negative lag and the observed 
spectrum. A somehow similar geometry (a standard 
disk with variable radius and a hotter, comptonizing 
cloud), although in a different context, has been suc- 
cessfully invoked to explain the spectral properties of 
Nova Muscae 1991 (e.g. Misra & Melia |1997j ; Esin, 
McClintock & Narayan [1997]). 
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In the next section we show that quantitative es- 
timates, based on this simple picture, are indeed in 
agreement with the observed colors and lags for rea- 
sonable values of the parameters. In our model the 
inner radius changes only by a factor ~ 3, so, for 
the sake of simplicity, in the following we neglect the 
small (~ 1.5, see eq. Q) change of T with rj n . 

3. NUMERICAL ESTIMATES 

The calculation of the emerging spectrum proceeds 
in two steps and follows the standard procedure dis- 
cussed by Sunyaev & Titarchuk (1980). First, we 



compute the time-dependent comptonized spectrum 
for a given seed photon distribution by solving the 
Kompaneets equation in a semi-infinite homogeneous 
medium and then use an escape probability algo- 
rithm to obtain the time evolution of the spectrum 
emerging from a medium of finite depth. We intro- 
duce a dimensionless time u = {n e UTc9)t and en- 
ergy x = hv/kT, where n e is the electron density, 
<jt the Thomson cross-section and 9 = kT/m e c 2 . 
In the following all lengths are expressed in units of 
GM/c 2 . The photon occupation number n(u,x) is 
the solution of the Kompaneets equation with an ap- 
propriate initial condition n(0,x). For rj n > r sc the 
initial photon occupation number is that of a mul- 
titemperature blackbody (calculated includin g GR 
and hardening effects, e.g. Novikov & Thorne 1973; 



Shimura & Takahara 1995| ) truncated at ri n whereas, 



for ri n < r sc , it corresponds to the superposition of a 
Bose-Einstein distribution at Tjj and a multitemper- 
ature blackbody truncated at r sc . 

Since in our model the long-term evolution corre- 
sponds to a variation of the inner edge of the disk, 
the Kompaneets equation has been solved numeri- 
cally for a sequence of values of rj n . Once n(u,x) has 
been evaluated, the time evolution of the escaping 
photons is given by 



N(u, x) = x 2 n(u, x)P(u, t) 



(2) 



where N is proportional to the count rate and P(u, r) 
is the escape probability distribution. The explicit 
form of P( u, t) depends on the geometry (Sunyaev & 
Titarchuk |l980|) . In our calculations we have used a 
weighted sum of the expressions for a spherical cloud 
with either a central or a diffused photon source. 
The total count is proportional to the convolution 



F(x) 



x 2 n(u, x)P(u, t) du . 



(3) 



We note that eq. (y) gives exactly the spectrum pro- 
duced by a stationary source of photons. Therefore it 
will be used to compute the hardness ratios on time 
intervals substantially larger than typical variability 



timescales in the continuum. On the other hand, if 
one is interested in following the time evolution of 
the spectrum over timescales ~ I/^qpo the relevant 
quantity to consider is N(u,x). 

To reproduce the observed time evolution of the 
color in a given energy band we have convolved 
N(u, x) with the interstellar absorption and with the 
detector response function 



C(E 1 



E 2 



E 2 ) oc / N(u,x)exp[-N H a(x)]A(x)dx 

J Ei 

where Njj is the column density. The lags in the 
continuum are computed comparing the time evolu- 
tion of the the soft and hard colors, C(2 — 5) and 
C(5 — 13), as obtained from eq. (||). The hardness 
ratios HR1 and HR2, defined as C(5 - 13)/C(2 - 5) 
and C(13 — 60)/C(2 — 5), are calculated again start- 
ing from eq. (0) but this time replacing N(u, x) with 
F{x). 

The calculations presented here have been obtained 
for Tjj = 15 keV, T\y = 1.5 keV and r sc = 20; the 
black hole mass and angular momentum are M = 
6M & and a = 0.1M. The model with rj n = 16 cor- 
responds to the usual picture in which a hot corona 
(T = Tjj) of moderate depth (r = tjj ~ 3) upscat- 
ters the soft photons produced by the disk; under 
these conditions the energy exchange in the warm 
region can be reasonably assumed to be negligible. 
The opposite situation is represented by the model 
with ri n = 6. Now the disk inner region lies entirely 
within the hot part of the corona which, because of 
the increase in the density, has tjj ^ 100. Escaping 
photons are then reprocessed in the outer, warm part 
of the corona (T = Tyy, r = Tw ~ 10). Intermedi- 
ate cases are more difficult to model in the present 
framework, because photons scatter in both regions 
and this would require a treatment of comptonization 
in a non-uniform medium. To avoid undue compli- 
cations, at this stage we assume that the effects of a 
temperature profile can be roughly described in terms 
of an "effective" temperature T, Ty/ < T < Tjj- For 
these models T and t were chosen in such a way to 
maintain y ~ 1. 

Our results are illustrated in Figure 0. As shown 
in Figure [|a, the observed hardness ratios are quite 
well reproduced by the model for a value of the col- 
umn density, Njj = 2 x 10 22 cm -2 , consistent with 



observations (Belloni et al. 2000 ). When the inner 
edge of the disk is far out, the spectrum is a power- 
law produced by thermal comptonization of soft disk 
photons. As ri n decreases, downscattering becomes 
progressively more important with respect to upscat- 
tering in the hot corona, the spectrum softens and, 
at the same time, the phase lag of the continuum de- 
creases. Eventually, systematic downscattering dom- 
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inates and soft (2-5 keV) photons, which underwent 
more scatterings, escape later than hard (5-13 keV) 
photons, causing an inversion in the sign of the lag. 
Figures ||b and ^c show the lag of the continuum as 
a function of HR1 and QPO frequency (the scatter- 
ing of the computed points is partially due to the 
poor sampling of the parameter space). Time lags 
have been computed assuming that the typical size 
of the comptonizing corona varies with time in ac- 
cordance with ri n in the range ~20-400. The rela- 
tion between vqpo and rj n was derived associating 
the model with rj n = 6 (16) to the observation with 
the largest (smallest) QPO frequency, 6 (0.6) Hz and 
yields 

(5) 



-0.43 
QPO ■ 



4. DISCUSSION 

In this Letter we have shown that a simple Comp- 
ton model can account for the main observational 
properties of GRS 1915+105 in state C. In partic- 
ular, the fundamental correlations (color-color, lag- 
color and lag-QPO frequency) are successfully repro- 
duced for a reasonable choice of the model parame- 
ters. Although we have adopted here a quite definite 
picture for the accretion flow, it is important to stress 
that the basic results are largely independent of the 
specific choice for the geometry. The decrease of the 
lag in the continuum and the softening of the spec- 
trum as the disk fills in are simply a consequence of 
the increase of the coronal depth and of the decrease 
of the temperature where radiation scatters before 
escaping. Thus, only three main ingredients are re- 
quired to account for the observed correlations within 
the framework of a simple thermal comptonization 
model: a non-isothermal corona of varying optical 
depth, a truncated standard disk and a source of 
more energetic (~ 10 keV) primary photons at small 
radii. Further assumptions (like the variable size of 
the corona and the existence of an "effective" comp- 
tonization temperature) are introduced only to match 
the observed values. It should be noted, however, 
that other geometrical settings for the time-varying 
accretion flow may work equally well. 

On the other hand, albeit not directly supported 
by hydrodynamics computations, the flow structure 
adopted here exhibits a certain degree of consistency 
in the way in which the QPO frequency and the coro- 



nal optical depth vary with inner disk radius. In 
fact, without entering into the mechanism that orig- 
inates the observed QPO, we have assumed that it 
is somehow produced at the inner edge of the disk 
and that its frequency varies with rj n . It is inter- 
esting to note that the relation between vqpo and 
ri n derived in § || (eq. Iff]) is similar to that found 
by Di Matteo & Psaltis ( |l999j ) combining the empir- 
ical relation of Psaltis, Belloni & van der Klis (1999) 
with the upper bound provided by the keplerian or- 
bital frequency. We stress, however, that our result 
has been obtained in an independent way and is a di- 
rect consequence of the values of the inner disk radius 
needed to match the observed colors and time lags. 
The contraction of the corona (a factor ~ 20) required 
to reproduce the largest /smallest lags implies an in- 
crease in the density (p oc i? -3 at constant coronal 
mass) consistent with the assumed ratio of the scat- 
tering depths, T H (r in = 6)/T H (r in = 16) ~ 100. The 
smallest value of the inner disk radius required by 
the model, rj n ~ 6M, is consistent with the assumed 
value a/M = 0.1. If the black hole is rapidly rotat- 
ing, as claimed by Zhang, Cui & Chen ( 1997 ) and 
Cui, Zhang & Chen (1998), then in low luminosity 
state C the disk does not reach the innermost stable 
circular orbit. 

A complete discussion of the behaviour of GRS 
1915+105 in states A and B is outside the scope of 
the present investigation. Here, we simply mention 
that models with high coronal optical depth but hard- 
photon starved (because, e.g., of the final collapse of 
the innermost, denser part of the corona) have hard- 
ness ratios in agreement with those observed in states 
A and B. Although no QPO has been observed in 
these states as yet, the positive detection of time lags 
correlated with colors in states A and B may lend 
further support to the present Compton model. 

Very recently, complex phase lag behavior in XTE 
J1550-564 has been reported by Cui, Zhang & Chen 
( pOOOl ). The power spectra are very similar to those 
in Reig et al. ( |2000| ) . However, the QPO phase-lag 
dependence is different. The fundamental/first har- 
monic shows a negative/zero lag, increasing in mag- 
nitude as the spectrum softens. In the light of the 
higher complexity of the timing properties of XTE 
J1550-564, any attempt to extend the present model 
to this source is premature. 
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Fig. 1. — A schematic view of the accretion flow in the central region of GRS 1915+105 in state C. The three panels illustrate 
the evolution as the disk rim moves closer to the hole (from left to right). Temperature increases from black to white. 
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Fig. 2. — a) color-color diagram for the model discussed in the text (filled diamonds) compared with observations (open circles); 
b) predicted and observed phase lags of the continuum versus HR1; c) predicted and observed phase lags of the continuum versus 
QPO frequency. Each model is labeled by the value of r^ n in units of GM/c 2 . 



